Genomic Insights into Mollusk Terrestrialization: Parallel and Convergent Gene Family Expansions as Key Facilitators in Out-of-the-Sea Transitions.

Animals abandoned their marine niche and successfully adapted to life on land multiple times throughout evolution, providing a rare opportunity to study the mechanisms driving large scale macroevolutionary convergence. However, the genomic factors underlying this process remain largely unknown. Here, we investigate the macroevolutionary dynamics of gene repertoire evolution during repeated transitions out of the sea in mollusks, a lineage that has transitioned to freshwater and terrestrial environments multiple independent times. Through phylogenomics and phylogenetic comparative methods, we examine ∼100 genomic data sets encompassing all major molluskan lineages. We introduce a conceptual framework for identifying and analyzing parallel and convergent evolution at the orthogroup level (groups of genes derived from a single ancestral gene in the species in question) and explore the extent of these mechanisms. Despite deep temporal divergences, we found that parallel expansions of ancient gene families played a major role in facilitating adaptation to nonmarine habitats, highlighting the relevance of the preexisting genomic toolkit in facilitating adaptation to new environments. The expanded functions primarily involve metabolic, osmoregulatory, and defense-related systems. We further found functionally convergent lineage-exclusive gene gains, while family contractions appear to be driven by neutral processes. Also, genomic innovations likely contributed to fuel independent habitat transitions. Overall, our study reveals that various mechanisms of gene repertoire evolution-parallelism, convergence, and innovation-can simultaneously contribute to major evolutionary transitions. Our results provide a genome-wide gene repertoire atlas of molluskan terrestrialization that paves the way toward further understanding the functional and evolutionary bases of this process.

Parallel and convergent genomic changes underlie independent subterranean colonization across beetles.

Adaptation to life in caves is often accompanied by dramatically convergent changes across distantly related taxa, epitomized by the loss or reduction of eyes and pigmentation. Nevertheless, the genomic underpinnings underlying cave-related phenotypes are largely unexplored from a macroevolutionary perspective. Here we investigate genome-wide gene evolutionary dynamics in three distantly related beetle tribes with at least six instances of independent colonization of subterranean habitats, inhabiting both aquatic and terrestrial underground systems. Our results indicate that remarkable gene repertoire changes mainly driven by gene family expansions occurred prior to underground colonization in the three tribes, suggesting that genomic exaptation may have facilitated a strict subterranean lifestyle parallelly across beetle lineages. The three tribes experienced both parallel and convergent changes in the evolutionary dynamics of their gene repertoires. These findings pave the way towards a deeper understanding of the evolution of the genomic toolkit in hypogean fauna.

Parallel duplication and loss of aquaporin-coding genes during the "out of the sea" transition as potential key drivers of animal terrestrialization.

One of the most important physiological challenges animals had to overcome during terrestrialization (i.e., the transition from sea to land) was water loss, which alters their osmotic and hydric homeostasis. Aquaporins are a superfamily of membrane water transporters heavily involved in osmoregulatory processes. Their diversity and evolutionary dynamics in most animal lineages remain unknown, hampering our understanding of their role in marine-terrestrial transitions. Here, we interrogated aquaporin gene repertoire evolution across the main terrestrial animal lineages. We annotated aquaporin-coding genes in genomic data from 458 species from seven animal phyla where terrestrialization episodes occurred. We then explored aquaporin gene evolutionary dynamics to assess differences between terrestrial and aquatic species through phylogenomics and phylogenetic comparative methods. Our results revealed parallel aquaporin-coding gene duplications during the ecological transition from marine to nonmarine environments (e.g., brackish, freshwater and terrestrial), rather than from aquatic to terrestrial ones, with some notable duplications in ancient lineages. In contrast, we also recovered a significantly lower number of superaquaporin genes in terrestrial arthropods, suggesting that more efficient oxygen homeostasis in land arthropods might be linked to a reduction in this type of aquaporin. Our results thus indicate that aquaporin-coding gene duplication and loss might have been one of the key steps towards the evolution of osmoregulation across animals, facilitating the "out of the sea" transition and ultimately the colonization of land.

Principled, practical, flexible, fast: a new approach to phylogenetic factor analysis.

Biological phenotypes are products of complex evolutionary processes in which selective forces influence multiple biological trait measurements in unknown ways. Phylogenetic comparative methods seek to disentangle these relationships across the evolutionary history of a group of organisms. Unfortunately, most existing methods fail to accommodate high-dimensional data with dozens or even thousands of observations per taxon. Phylogenetic factor analysis offers a solution to the challenge of dimensionality. However, scientists seeking to employ this modeling framework confront numerous modeling and implementation decisions, the details of which pose computational and replicability challenges.We develop new inference techniques that increase both the computational efficiency and modeling flexibility of phylogenetic factor analysis. To facilitate adoption of these new methods, we present a practical analysis plan that guides researchers through the web of complex modeling decisions. We codify this analysis plan in an automated pipeline that distills the potentially overwhelming array of decisions into a small handful of (typically binary) choices.We demonstrate the utility of these methods and analysis plan in four real-world problems of varying scales. Specifically, we study floral phenotype and pollination in columbines, domestication in industrial yeast, life history in mammals, and brain morphology in New World monkeys.General and impactful community employment of these methods requires a data scientific analysis plan that balances flexibility, speed and ease of use, while minimizing model and algorithm tuning. Even in the presence of non-trivial phylogenetic model constraints, we show that one may analytically address latent factor uncertainty in a way that (a) aids model flexibility, (b) accelerates computation (by as much as 500-fold) and (c) decreases required tuning. These efforts coalesce to create an accessible Bayesian approach to high-dimensional phylogenetic comparative methods on large trees.

Endocranial asymmetry in New World monkeys: a comparative phylogenetic analysis of morphometric data.

Brain lateralization is a widespread phenomenon although its expression across primates is still controversial due to the reduced number of species analyzed and the disparity of methods used. To gain insight into the diversification of neuroanatomical asymmetries in non-human primates we analyze the endocasts, as a proxy of external brain morphology, of a large sample of New World monkeys and test the effect of brain size, home range and group sizes in the pattern and magnitude of shape asymmetry. Digital endocasts from 26 species were obtained from MicroCT scans and a set of 3D coordinates was digitized on endocast surfaces. Results indicate that Ateles, Brachyteles, Callicebus and Cacajao tend to have a rightward frontal and a leftward occipital lobe asymmetry, whereas Aotus, Callitrichinae and Cebinae have either the opposite pattern or no directional asymmetry. Such differences in the pattern of asymmetry were associated with group and home range sizes. Conversely, its magnitude was significantly associated with brain size, with larger-brained species showing higher inter-hemispheric differences. These findings support the hypothesis that reduction in inter-hemispheric connectivity in larger brains favors the lateralization and increases the structural asymmetries, whereas the patterns of shape asymmetry might be driven by socio-ecological differences among species.

Characterizing and Comparing Phylogenetic Trait Data from Their Normalized Laplacian Spectrum.

The dissection of the mode and tempo of phenotypic evolution is integral to our understanding of global biodiversity. Our ability to infer patterns of phenotypes across phylogenetic clades is essential to how we infer the macroevolutionary processes governing those patterns. Many methods are already available for fitting models of phenotypic evolution to data. However, there is currently no comprehensive nonparametric framework for characterizing and comparing patterns of phenotypic evolution. Here, we build on a recently introduced approach for using the phylogenetic spectral density profile (SDP) to compare and characterize patterns of phylogenetic diversification, in order to provide a framework for nonparametric analysis of phylogenetic trait data. We show how to construct the SDP of trait data on a phylogenetic tree from the normalized graph Laplacian. We demonstrate on simulated data the utility of the SDP to successfully cluster phylogenetic trait data into meaningful groups and to characterize the phenotypic patterning within those groups. We furthermore demonstrate how the SDP is a powerful tool for visualizing phenotypic space across traits and for assessing whether distinct trait evolution models are distinguishable on a given empirical phylogeny. We illustrate the approach in two empirical data sets: a comprehensive data set of traits involved in song, plumage, and resource-use in tanagers, and a high-dimensional data set of endocranial landmarks in New World monkeys. Considering the proliferation of morphometric and molecular data collected across the tree of life, we expect this approach will benefit big data analyses requiring a comprehensive and intuitive framework.

Understanding the effect of competition during evolutionary radiations: an integrated model of phenotypic and species diversification.

Competition can drive macroevolutionary change, for example during adaptive radiations. However, we still lack a clear understanding of how it shapes diversification processes and patterns. To better understand the macroevolutionary consequences of competition, as well as the signal left on phylogenetic data, we developed a model linking trait evolution and species diversification in an ecological context. We find four main results: first, competition spurs trait diversity but not necessarily species richness; second, competition produces slowdowns in species diversification even in the absence of explicit ecological limits, but not in phenotypic diversification even in the presence of such limits; third, early burst patterns do not provide a reliable way of testing for adaptive radiations; and fourth, looking for phylogenetic signal in trait data and support for phenotypic models incorporating competition is a better alternative. Our results clarify the macroevolutionary consequences of competition and could help design more powerful tests of adaptive radiations in nature.

Cranial and endocranial diversity in extant and fossil atelids (Platyrrhini: Atelidae): A geometric morphometric study.

OBJECTIVES: Platyrrhines constitute a diverse clade, with the modern Atelidae exhibiting the most variation in cranial and endocast morphology. The processes responsible for this diversification are not well understood. Here, we present a geometric morphometric study describing variation in cranial and endocranial shape of 14 species of Alouatta, Ateles, Brachyteles, and Lagothrix and two extinct taxa, Cartelles and Caipora. METHODS: We examined cranial and endocranial shape variation among species using images reconstructed from CT scans and geometric morphometric techniques based on three-dimensional landmarks and semilandmarks. Principal components analyses were used to explore variation, including the Procrustes shape coordinates, summing the logarithm of the Centroid Size, the common allometric component, and residual shape components. RESULTS: Differences in endocranial shape are related to a relative increase or decrease in the volume of the neocortex region with respect to brainstem and cerebellum regions. The relative position of the brainstem varies from a posterior position in Alouatta to a more ventral position in Ateles. The shape of both the cranium and endocast of Caipora is within the observed variation of Brachyteles. Cartelles occupies the most differentiated position relative to the extant taxa, especially in regards to its endocranial shape. CONCLUSIONS: The pattern of variation in the extant species in endocranial shape is similar to the variation observed in previous cranial studies, with Alouatta as an outlier. The similarities between Caipora and Brachyteles were unexpected and intriguing given the frugivorous adaptations inferred from the fossil's dentition. Our study shows the importance of considering both extant and fossil species when studying diversification of complex traits.

A Penalized Likelihood Framework for High-Dimensional Phylogenetic Comparative Methods and an Application to New-World Monkeys Brain Evolution.

Working with high-dimensional phylogenetic comparative data sets is challenging because likelihood-based multivariate methods suffer from low statistical performances as the number of traits $p $ approaches the number of species $n $ and because some computational complications occur when $p $ exceeds $n$. Alternative phylogenetic comparative methods have recently been proposed to deal with the large $p $ small $n $ scenario but their use and performances are limited. Herein, we develop a penalized likelihood (PL) framework to deal with high-dimensional comparative data sets. We propose various penalizations and methods for selecting the intensity of the penalties. We apply this general framework to the estimation of parameters (the evolutionary trait covariance matrix and parameters of the evolutionary model) and model comparison for the high-dimensional multivariate Brownian motion (BM), Early-burst (EB), Ornstein-Uhlenbeck (OU), and Pagel's lambda models. We show using simulations that our PL approach dramatically improves the estimation of evolutionary trait covariance matrices and model parameters when $p$ approaches $n$, and allows for their accurate estimation when $p$ equals or exceeds $n$. In addition, we show that PL models can be efficiently compared using generalized information criterion (GIC). We implement these methods, as well as the related estimation of ancestral states and the computation of phylogenetic principal component analysis in the R package RPANDA and mvMORPH. Finally, we illustrate the utility of the new proposed framework by evaluating evolutionary models fit, analyzing integration patterns, and reconstructing evolutionary trajectories for a high-dimensional 3D data set of brain shape in the New World monkeys. We find a clear support for an EB model suggesting an early diversification of brain morphology during the ecological radiation of the clade. PL offers an efficient way to deal with high-dimensional multivariate comparative data.

Multiple factors behind early diversification of skull morphology in the continental radiation of New World monkeys.

Understanding the origin of diversity is a fundamental problem in evolutionary biology. The null expectation for the evolutionary diversification is that all changes in biological diversity are the result of random processes. Adaptive radiations depart from this expectation as ecological factors and natural selection are supposed to play a central role in driving exceptional diversification. However, it is not well understood how large-scale continental radiations, given their characteristics, fit to these opposing theoretical models. Here, we used phylogenetic comparative methods and geometric morphometrics to study the evolutionary process of cranial diversification in the continental radiation of New World monkeys. Particularly, we tested several alternative evolutionary scenarios for morphological evolution in the clade. Results indicated that despite the platyrrhine radiation being old and geographically widespread, the formative patterns arising from the initial stages of diversification probably associated with an adaptive radiation can still be recognized today. We also show that no single explored factor (e.g., ecological or allometric) can be invoked as a complete explanation for the observed phenotypic diversity patterns in the clade and, moreover, that different cranial regions exhibit particular macroevolutionary patterns. Together, our results highlight the evident complexity behind large-scale evolutionary radiations.

Early evolutionary diversification of mandible morphology in the New World monkeys (Primate, Platyrrhini).

New World monkeys (order Primates) are an example of a major mammalian evolutionary radiation in the Americas, with a contentious fossil record. There is evidence of an early platyrrhine occupation of this continent by the Eocene-Oligocene transition, evolving in isolation from the Old World primates from then on, and developing extensive morphological and size variation. Previous studies postulated that the platyrrhine clade arose as a local version of the Simpsonian ecospace model, with an early phase involving a rapid increase in morphological and ecological diversity driven by selection and ecological opportunity, followed by a diversification rate that slowed due to niche-filling. Under this model, variation in extant platyrrhines, in particular anatomical complexes, may resemble patterns seen among middle-late Miocene (10-14 Ma) platyrrhines as a result of evolutionary stasis. Here we examine the mandible in this regard, which may be informative about the dietary and phylogenetic history of the New World monkeys. Specifically, we test the hypothesis that the Simpsonian ecospace model applies to the platyrrhine mandible through a geometric morphometric analysis of digital images of the jaws of extant and extinct species, and we compare these results to those obtained using a phylogenetic comparative approach based on extant species. The results show a marked phylogenetic structure in the mandibular morphology of platyrrhines. Principal component analyses highlight the morphological diversity among modern forms, and reveal a similar range of variation for the clade when fossil specimens are included. Disparity-Through-Time analysis shows that most of the shape variation between platyrrhines originated early in their evolution (between 20 and 15 Ma). Our results converge with previous studies of body mass, cranial shape, the brain and the basicranium to show that platyrrhine evolution might have been shaped by an early increase in morphological variation followed by a decelerated rate of diversification and evolutionary stasis.

Brain shape convergence in the adaptive radiation of New World monkeys.

Primates constitute one of the most diverse mammalian clades, and a notable feature of their diversification is the evolution of brain morphology. However, the evolutionary processes and ecological factors behind these changes are largely unknown. In this work, we investigate brain shape diversification of New World monkeys during their adaptive radiation in relation to different ecological dimensions. Our results reveal that brain diversification in this clade can be explained by invoking a model of adaptive peak shifts to unique and shared optima, defined by a multidimensional ecological niche hypothesis. Particularly, we show that the evolution of convergent brain phenotypes may be related to ecological factors associated with group size (e.g., social complexity). Together, our results highlight the complexity of brain evolution and the ecological significance of brain shape changes during the evolutionary diversification of a primate clade.

Encephalization and diversification of the cranial base in platyrrhine primates.

The cranial base, composed of the midline and lateral basicranium, is a structurally important region of the skull associated with several key traits, which has been extensively studied in anthropology and primatology. In particular, most studies have focused on the association between midline cranial base flexion and relative brain size, or encephalization. However, variation in lateral basicranial morphology has been studied less thoroughly. Platyrrhines are a group of primates that experienced a major evolutionary radiation accompanied by extensive morphological diversification in Central and South America over a large temporal scale. Previous studies have also suggested that they underwent several evolutionarily independent processes of encephalization. Given these characteristics, platyrrhines present an excellent opportunity to study, on a large phylogenetic scale, the morphological correlates of primate diversification in brain size. In this study we explore the pattern of variation in basicranial morphology and its relationship with phylogenetic branching and with encephalization in platyrrhines. We quantify variation in the 3D shape of the midline and lateral basicranium and endocranial volumes in a large sample of platyrrhine species, employing high-resolution CT-scans and geometric morphometric techniques. We investigate the relationship between basicranial shape and encephalization using phylogenetic regression methods and calculate a measure of phylogenetic signal in the datasets. The results showed that phylogenetic structure is the most important dimension for understanding platyrrhine cranial base diversification; only Aotus species do not show concordance with our molecular phylogeny. Encephalization was only correlated with midline basicranial flexion, and species that exhibit convergence in their relative brain size do not display convergence in lateral basicranial shape. The evolution of basicranial variation in primates is probably more complex than previously believed, and understanding it will require further studies exploring the complex interactions between encephalization, brain shape, cranial base morphology, and ecological dimensions acting along the species divergence process.

Modeling lineage and phenotypic diversification in the New World monkey (Platyrrhini, Primates) radiation.

Adaptive radiations that have taken place in the distant past can now be more thoroughly studied with the availability of large molecular phylogenies and comparative data drawn from extant and fossil species. Platyrrhines are a good example of a major mammalian evolutionary radiation confined to a single continent, involving a relatively large temporal scale and documented by a relatively small but informative fossil record. Here, we present comparative evidence using data on extant and fossil species to explore alternative evolutionary models in an effort to better understand the process of platyrrhine lineage and phenotypic diversification. Specifically, we compare the likelihood of null models of lineage and phenotypic diversification versus various models of adaptive evolution. Moreover, we statistically explore the main ecological dimension behind the platyrrhine diversification. Contrary to the previous proposals, our study did not find evidence of a rapid lineage accumulation in the phylogenetic tree of extant platyrrhine species. However, the fossil-based diversity curve seems to show a slowdown in diversification rates toward present times. This also suggests an early high rate of extinction among lineages within crown Platyrrhini. Finally, our analyses support the hypothesis that the platyrrhine phenotypic diversification appears to be characterized by an early and profound differentiation in body size related to a multidimensional niche model, followed by little subsequent change (i.e., stasis).

Divergence Times and the Evolutionary Radiation of New World Monkeys (Platyrrhini, Primates): An Analysis of Fossil and Molecular Data.

The estimation of phylogenetic relationships and divergence times among a group of organisms is a fundamental first step toward understanding its biological diversification. The time of the most recent or last common ancestor (LCA) of extant platyrrhines is one of the most controversial among scholars of primate evolution. Here we use two molecular based approaches to date the initial divergence of the platyrrhine clade, Bayesian estimations under a relaxed-clock model and substitution rate plus generation time and body size, employing the fossil record and genome datasets. We also explore the robustness of our estimations with respect to changes in topology, fossil constraints and substitution rate, and discuss the implications of our findings for understanding the platyrrhine radiation. Our results suggest that fossil constraints, topology and substitution rate have an important influence on our divergence time estimates. Bayesian estimates using conservative but realistic fossil constraints suggest that the LCA of extant platyrrhines existed at ca. 29 Ma, with the 95% confidence limit for the node ranging from 27-31 Ma. The LCA of extant platyrrhine monkeys based on substitution rate corrected by generation time and body size was established between 21-29 Ma. The estimates based on the two approaches used in this study recalibrate the ages of the major platyrrhine clades and corroborate the hypothesis that they constitute very old lineages. These results can help reconcile several controversial points concerning the affinities of key early Miocene fossils that have arisen among paleontologists and molecular systematists. However, they cannot resolve the controversy of whether these fossil species truly belong to the extant lineages or to a stem platyrrhine clade. That question can only be resolved by morphology. Finally, we show that the use of different approaches and well supported fossil information gives a more robust divergence time estimate of a clade.